arXiv:hep-ph/9808206vl 3 Aug 1998 


IHEP 98-41 
.June 1998 


THE MODEL FOR QCD RUNNING COUPLING 
CONSTANT WITH DYNAMICALLY GENERATED 
MASS AND ENHANCEMENT IN THE INFRARED 

REGIONf] 


Aleksey I. Alekseev 

Institute for High Energy Physics, Protvino, Moscow Region, 142284 Russia 

E-mail: alekseev@mx.ihep.su 


Abstract 

Nonperturbative studies of the strong running coupling constant 
in the infrared region are discussed. Starting from the analyses of 
the Dyson - Schwinger equations in the gauge sector of QCD, the 
conclusion is made on an incomplete fixing of the perturbation the¬ 
ory summation ambiguity within ’’(forced) analytization procedure” 

(called also a dispersive approach). A minimal model for a s (q 1 2 ) is 
proposed so that the perturbative time-like discontinuity is preserved 
and nonperturbative terms not only remove the Landau singularity 
but also provide the ultraviolet convergence of the gluon condensate. 
Within this model, on the one hand, the gluon zero modes are en¬ 
hanced (the dual superconductor property of the QCD vacuum) and, 
on the other hand, dynamical gluon mass generation is realized, with 
m g estimated as 0 .6GeV. The uncertainty connected with the division 
into perturbative and nonperturbative contributions is discussed with 
the gluon condensate taken as an example. 

1 Talk presented at the Workshop on Methods in Non-Perturbative Field Theory, Febru¬ 

ary 2 - February 13, 1998, Adelaide, Australia 



1 Introduction 


The infrared region corresponding to large distances is of exceptional interest 
by virtue of its responsibility for the confinement mechanism and its inac¬ 
cessibility to perturbative methods. A great number of papers are devoted 
to the infrared behaviour of gluon Green functions but a commonly accepted 
opinion is absent in the literature. 

Recently a possibility of power infrared behaviour of the gluon propagator 
in the covariant Landau gauge has been discussed, B0 with a partial account 
for ghosts. It was found that in the approximation considered, the gluon 
propagator at small values of Euclidean momenta vanished (in this case, 
according to gauge identities, all the gluon vertexes turn out to be singular), 
the ghost propagator was singular, and a running coupling constant had a 
large but finite value at zero. 

In the axial gauge in the framework of nonperturbative Baker - Ball - 
Zachariasen (BBZ) approach, || which seems to be adequate to discuss the 
possibility of the infrared behaviour not being too singular, the problem of 
consistency of the behaviour D(q) ~ ( q 2 )~ c , q 2 —> 0 was considered. It has 
been shown [(]]. |5|] that for quite a wide interval of non-integer (non-half- 
integer) values of c, — 1 < c < 3, there are no solutions. The possibility of 
’’frozen” interaction in the infrared region was also considered in the frame¬ 
work of the above approach and the result was negative as well. || 

On the other hand, the enhanced infrared behaviour of the gluon propa¬ 
gator, of the form D(q) ~ l/(g 2 ) 2 , q 2 —> 0 is physically motivated, useful in 
applications (e.g. in the quark sector), and it was obtained in a number of 
different approaches. 0, [| This behaviour asymptotically solves the Dyson 
- Schwinger equation in the axial gauge, but with this we should decline 
one of the basic assumptions of the BBZ approach and take into account a 
transverse part of the triple gluon vertex. || Accordingly, we should decline 
an appealing iteration scheme [pj to find solutions of the Dyson - Schwinger 
equations for the higher Green functions. 

Versions of the infrared behaviour of the gluon propagator mentioned 
above do not exhaust all the possibilities that can be realized beyond per¬ 
turbation theory. 
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2 Model for the Running Coupling Constant 
with Nonperturbative Contributions Sup¬ 
pressed in the Ultraviolet Region 


Let us consider the model for a s introduced in Refs. [[| and discussed fur¬ 
ther in Ref. Q We begin with ’’analytized” expression for the QCD running 
coupling constant [jTOj 

A 2 


a 


(i) 


/ 2n 4?r 

(0 = 7 - 
Oo 


+ 


ln(g 2 /A 2 ) A 2 — q 


( 1 ) 


obtained with the use of the idea [|ll|] on the cancellation of the ghost singular¬ 
ity by nonperturbative contributions. We see that this expression has a non¬ 
perturbative tail, with the behaviour 1/q 2 at q 2 —> oo. To answer the question 
whether this behaviour is admissible, let us consider an important physical 
quantity, namely, the gluon condensate K = < vac \ a s /iT : F“ u F° u :| vac > . 
Up to the quadratic approximation in the gluon fields, the gluon condensate 
K is defined by nonperturbative contributions in the transverse part of the 
gluon propagator. Normal ordering of the operators product is defined with 
respect to perturbative vacuum |pert >, as the averaging in the expression 
for the condensate is carried out with true physical vacuum state |vac >. 
The following chain of equations describes the method by which we obtain 
closed expression for the gluon condensate: 


QL 

K = Jim < vac\-^ : F“„(x)F“„(y) : \vac (2) 

« 2 8 ab lim < uac|— T[(<9 M A“(x) — d v A a (x))d^A b v (y)\\vac > —pert = (3) 

*y tt 

= 2 5 ab — lim (5^d x x d y x - d x u d v )N~ l f dAA a (x)A b v (y)e s[A]+JA | J=0 -pert = 
7T x ~*y J 

(4) 

cy 

= 2lim (S,„dld\ - {D%(x, y) - D^(x, y)]= (5) 

7T x ^y 

= 2^1im(VS5S?-S;^)/ ■^Y^ ,k{x - y) [D^{k)-D^(k)] = (6) 
= 2 J ^L(fc 2 V - k.K)\D„„(k) - D%T(k)] = (7) 
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= jdk(k%„ - k„K)D^(k)[Z(k) - Z^ H (k)} = ( 8 ) 

= ^-Jdk[Z(k)-Z-^(k)} = (9) 

= 4 / ik\a,(k?) - af H (k 2 )] = 4 ; jiyyd*°"'‘ n {y). ( 10 ) 


We shall use Eq. (|T0|) to evaluate the gluon condensate (§). The transi¬ 
tion from Eq. (17]) to (|8|) is made with the assumption of proportionality of 
the nonperturbative part of the propagator to the free one. The projector in 
Eq. (^) removes the terms of the free axial gauge propagator D^{k) which 
depend on the gauge vector. Assuming Z is independent of gauge parameter 
y = (k7]) 2 /k 2 7] 2 : the employment of the axial gauge lets one obtain Eq. (0) 
from (||), thereby expressing the gluon condensate in terms of nonperturba¬ 
tive part of running coupling constant a s (k 2 ). The one-loop ’’analytized” 
behaviour of Eq. ([[]) leads to a quadratic divergence in Eq. flip] ) at infinity 
and this is true for the two- and three-loop expressions |L(| of the analytiza- 
tion approach. According to the results of Refs. 0 ,&§ let us add in Eq. ( 0 ) 
the isolated infrared singular term of the form 1/q 2 . This term is harmless at 
zero and it can improve the behaviour of the integrand at infinity and make 
the integral logarithmical divergent. To make the integral (|10]) convergent at 
infinity, it is sufficient to add one more isolated singular term of a pole type 
with parameters chosen appropriately. In this sense the model we come to 
is minimal. The expression, we obtain for the running coupling constant, is 
the following: 


47T / 1 A 2 cA 2 (1 —c)A 2 \ 

b 0 \ln(q 2 / A 2 ) A 2 — q 2 q 2 q 2 + m 2 J ’ 


( 11 ) 


with mass parameter m g = A j\Jc — 1. It is worth noting that an account of 
nonperturbative contributions in Eq. ( ]TT| ) preserves a perturbative time-like 
cut of Eq. (]T|). With the given value of the QCD scale parameter A, the 
parameter c can be fixed by the string tension k, assuming the linear con¬ 
finement V(r) ~ nr = a 2 r at r —> oo. We define the potential V(r) of static 
qq interaction by means of three-dimensional Fourier transform of Oi s (q 2 )/q 2 
with the contributions of only one dressed gluon exchange taken into account. 
This gives the relation cA 2 = (36o/87r)a 2 . Taking a ~ 0.42 GeV, one obtains 
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c = A 2 /A 2 where Aj = 3b 0 n/8n ~ (0.434 GeV) 2 (b 0 = 9 in the case of 3 light 
flavours). For m g one obtains 

m, = A 2 / A 2 - 73 , (12) 

and the tachion absence condition limits the parameter A to A < 434 MeV. 

ft is seen from Eq. that the pole singularities are situated at two 
points q 2 = 0 and q 2 = —m 2 . This corresponds to the two effective gluon 
masses, 0 and rn g . Therefore, the physical meaning of the parameter rn g is 
not the constituent gluon mass, but rather the mass of the exited state of 
the gluon. 


3 Gluon Condensate and Nonperturbative Scale 


The acceptance of the cancellation mechanism for the nonphysical pertur¬ 
bation theory singularity by the nonperturbative contributions leads to the 
necessity of a supplementary definition of integral (|l0j) near point k 2 = A 2 . 
This problem can be reformulated as a problem of dividing the perturba¬ 
tive and nonperturbative contributions in a s resulting in the introduction of 
a parameter ko ~ 1 GeV. The following procedure is our definition of the 
regularized perturbative and nonperturbative parts of a s {k 2 )\ 

A 2 

+ 


ot s (k 2 ) = — 


where 



ln(fc 2 /A 2 ) A 2 -k 2 


cA 2 (1 — c)A 2 ' 

+ -tt- + V ' 


k 2 


k 2 + m 2 


(Fj+ar^'t* 2 ) = a?;’(k 2 ) + <C3jr‘(k 


(V) = 


47T / 


+ 


A 2 


b 0 \ln(k 2 /A 2 ) A 2 -k 2 

has no power corrections at k 2 —> oo and 
An ( A 2 


9{kl - k 2 


a 


nonpert(1^2} _ 

s reg 


(k 2 ) = 


b 0 \A 2 - k 2 


0(k 2 - k 2 ) + — + 


cA 2 (1 - c) A 2 


k 2 + m 2 


(13) 

(14) 

(15) 

(16) 


at k 2 —>■ oo has the same power corrections as o i ™ onpert (k 2 ), namely ~ A 6 /k 6 . 
Thus, we regularize ^ the usual perturbation theory in the infrared region. 

2 See also Ref. |I2| where the problem of perturbative and nonperturbative contributions 

to a s is discussed and the definition of infrared finite regularized perturbative part of a s 

is suggested. 
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Figure 1: The dependence of the gluon condensate fourth root Ji 1 / 4 on A. 
The parameter k 0 varies from 0.5 GeV (the lowest curve) up to 3.0 Geld (the 
highest one) at an interval of 0.25 Geld. The dashed line corresponds to the 
exceptional value k 0 = 0.777GeV. The ’’standard” level K 1 ! A = 0.33 Geld is 
indicated. 


The arbitrariness of the procedure is parametrized by a scale parameter ko- 
With account for the stated above, let us calculate the gluon condensate (p~Q|) . 
We obtain 


K(A,k 0 ) 




(17) 


In Fig. 1 the dependence of AT 1 / 4 on A is shown for different values of ko in 
the interval (0.5 — 3.0) GeV. For our estimates we shall use the ’’standard” 
value ]TB[ of the gluon condensate (0.33 GeV) 4 . We can see that at ko < ko — 
0.777 GeV gluon condensate (p~7|) is smaller than its standard value. It means 
that there exists a lower limit for the value of parameter ko . The value 
ko = ko turns out to be exceptional because with this choice only one value 
A = A ~ 375 MeV provides a necessary value of gluon condensate. In this 
case, according to ©. m g = rri g ~ 0.6 GeV. At k 0 > k 0 the two values of 
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A give the needed value of the gluon condensate. When k 0 increases one of 
them increases and tends to Ai and the other value of A decreases and, at 
k 0 > 2 GeV, becomes less than 150 MeV. 


4 Gluon Condensate and Vacuum Energy Den¬ 
sity 


Let us represent expression ([II]) in an explicitly renormalization invariant 
form, ft can be done without solving the differential renormalization group 
equations. In this order we write a s (q 2 ) = g 2 (q 2 /g 2 , g)/An and use the 
normalization condition g 2 (l, g) = g 2 . Then we obtain the equation for the 
wanted dependence of the parameter A 2 on g 2 and g 2 : 

1 A 2 A 2 (1 — c)A 2 " 

ln(p 2 /A 2 ) A 2 — /i 2 C g 2 g 2 + m 2 

For dimensional reasons A 2 = g 2 exp{ — ip(x)}, where x = bog 2 /lQn 2 = 
b 0 a s /4:7T, and for the function tp(x) we obtain the equation: 


2 4vr 
g 4vr = — 

On 


X = 


1 

<p(x) 




+ ce~^ x) - 


(o- l ) 2 

(c- l)e^G) + 1' 


The solution of this equation at c > 1 is function <p(x), which has the be¬ 
haviour <p(x) — 1/x at x —> 0 and <p(x) — —ln(x/c) at x —> +oo. The 
relation obtained ensures the renormalization invariance of a s (q 2 ). At low 
g 2 , we obtain A 2 = p 2 exp{—47 t/(6oQ: s )}, which indicates the essentially non- 
perturbative character of three last terms in Eq. (0 and these terms are 
absent in the usual perturbation theory. We will need further the /3-function 
which can be found by the equation 

= p(f) . (18) 

Here u = q 2 /g 2 ■ Differentiating the running coupling ([IT]) in u and assuming 
u = 1, we obtain 




167T 2 f 1 1 1 (c — l) 2 1 

bo \ <p(x) <p 2 {x) (e^G) — l ) 2 ((c — l)eV 3 ( a: ) -)- l) 2 J 

(19) 
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Figure 2: The [3 - function for c = 1, 2, 3,4, 5. 


which is simplified using equation for </?(x). Then, knowing the behaviour of 
<p(x) at x —> 0, oo, we can find 


/ Kg 2 ) * 


167T 2 


0 4 . 2 

9 + 9 


0, 


( 20 ) 


M) ^ -9 2 


167T 2 


c(c-2) + 0(l/lns 2 ), j 2 -» oo. 


( 21 ) 


Fig. 2 illustrates the dependence /3(g 2 ) for c = 1,2, ...,5. For all g 2 > 0 the 
/3-function is negative definite. Let us consider the trace anomaly [14] for the 
energy - momentum tensor of the gluon field 


B(q 2 ) 

0 = ■ F a F a • (22) 

yj im 2g 2 ' ^ ^ ' ' y 

According to the definition the vacuum is relativistically invariant. So, for 
the vacuum energy density we have e vac = (1/4) < vac\®^{vac >. For 
(3(g 2 ) = — 6 0 g 4 /(167r 2 ), b 0 = 9 and K— (0.33 GeV) 4 , the vacuum energy 
density is e vac ~ —(240 MeV) A . If we introduce quarks, we destroy the 
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Figure 3: K{fx)/K{fi —*• oo) = — (b 0 /16n 2 )g 4 //3(g 2 ) and (K (//)/ K(/j, —> 

oo)) 1 / 4 . 

nonperturbative vacuum in some region (bag). The vacuum inside the bag 
is perturbative and its energy density is zero. The difference of the vacuum 
energy densities inside and outside the bag is the reason of external pressure 
on the bag. 

It is important for us that because of /3 < 0, the maximum of gluon 
condensate corresponds to the minimum of nonperturbative vacuum energy 
density. For our model a s (q 2 ) it means that the values 

A = A ~ 375 MeV, ko = ko ~ 777 MeV, rn g = fh g ~ GOOMeV" (23) 

are advantageous from the viewpoint of the energy argument. Let us clarify 
what the dependence of the gluon condensate on renormalization parameter 
/i is. As a physical quantity, the vacuum energy density is independent of /j. 
The renormalized coupling constant g depends on /j,, g 2 / 47T = g 2 {l, g)/ = 
a s (g, 2 ), so as /i —» oo g 2 — » 0 and at /r —> 0 g 2 —> oo. Using (|20|) , (|2l|) , we 
can write 

W = < vac \~ : : \vac >= (24) 
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= < vac \^ : F* V F* V : |vac >| M ^oo= 

(25) 


= < vac \ ■ F lv F lv ■ \ vac >U-0 • 

(26) 

From (pi]) 

, (^) at bo = 9, we have the ratio 



K{n) 9 g 3 4 5 

K(n—> oo) 167T 2 /3(g 2 ) ’ 

(27) 


which describes the dependence of gluon condensate on /i. For the /3-function (|T9| ) 
in Fig. 3 by a solid line we show the ratio K(n)/K(ii —> cx)) as function of 
(i. The parameter c = 1, 3476 corresponds to the choice Ai = 434 MeV, 

A = A = 375 MeV. We can see that for /i > 1 GeV ratio (^7|) is practically 
unity. In the same figure the dependence > oo)) 1 / 4 on /i is 

shown by dots. 
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